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Abstract 


Quantifying the spacecraft charging risks and corresponding hazards for 
the International Space Station (ISS) requires a plasma environment 
specification describing the natural variability of ionospheric temperature 
and density. Empirical ionospheric specification and forecast models such as 
the International Reference Ionosphere (IRI) model typically only provide 
estimates of the long term (seasonal) mean electron temperature and density 
values for the low Earth orbit environment. Knowledge of the electron 
temperature and density variability as well as the likelihood of extreme 
deviations from the mean values are required to estimate both the magnitude 
and frequency of occurrence of potentially hazardous spacecraft charging 
environments for a given ISS construction stage and flight configuration. 

This paper describes the statistical analysis of historical ionospheric low Earth 
orbit plasma measurements used to estimate the variability of density (Ne) and 
temperature (Te) in the ISS flight environment. The statistical variability 
analysis of Ne and Te enables calculation of the expected frequency of 
occurrence of any particular values of Ne and Te, especially those 
that correspond to possibly hazardous spacecraft charging 
environments. The database used in the original analysis included 
measurements from the AE-C, AE-D, and DE-2 satellites. Recent work on the 
database has added additional satellites to the database and ground based 
incoherent scatter radar observations as well. Deviations of the data values 
from the International Reference Ionosphere estimated Ne, Te parameters for 
each data point provide a statistical basis for modeling the deviations of the 
plasma environment from the IRI model output. This technique, while 
developed specifically for the Space Station, can also be generalized to 


provide ionospheric plasma environment risk specification models for low 
Earth orbit over an altitude range of 200 km through approximately 1000 km. 
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Overview 


• Brief review of work to date 

• Current work 

- New data 

- Documentation 

- Forward word (schedule) 


Brief review of work to date 


Variability Study 




ISS Program requires estimate of Te, Ne variability: 

- range of possible values at ISS altitudes 

- statistical variations in parameters 

- values for sunrise (eclipse exit) at high latitude particularly 
important 

Ionosphere models (e.g., IRI-90, IRI-95) don’t contain information on 
variability at short time scales 

- IRI 2001 does contain modules to include geomagnetic storms 
Models are too “smooth” 

- data is required to obtain variability in environment about mean values 

Analysis of historical spacecraft data can provide estimate of 
variability (na or other statistical figure of merit) in ionospheric 
parameters 



4 



Data 


Spacecraft 

• AE-C 

- Circular, 68.1 deg inc 

Elliptical 150 km x 4300 km x 68.1 deg inc 
1973-12-16-00-00 - 1978-12-11-59-59 

• AE-D 

- Circular, 90. 1 deg inc 

- Elliptical, 154 km x 3816 km x 90.1 deg inc 

- 1975-10-06-00-00 - 1976-01-29-59-59 

• AE-E 

- Circular, 19.7 deg inc 

- Elliptical, 156 km x 2983 km x 90.1 deg inc 

- 1975-12-01-00-00 - 1981-06-06 

• DE 2 

- Elliptical, 300 km x 1000 km x 89.99 deg 
1981-08-06-00-00 - 1983-02-15-59-59 

• Explorer 22 

- 889 km x 1081 km x 79.7 deg 

- 1964-10-10 to 1967-08-26 

• Explorer 31 

- 500 km x 3000 km x 79.8 deg 

- 1965-10-29 to 1968-08-20 


Data Used in Study 

[1973-12-16 to 1978-12-11] 

608,139 records* 


[1975-10-06 to 1976-01-29] 

100,395 records* 


[1975-12-01 to 1981-05-24] 

464,690 records* 


[1981-08-06 to 1983-02-15] 

693,112 records* 

[1964-10-10 to 1967-08-26] 

37,476 records* 

[1965-10-29 to 1968-08-20] 

55,870 records* 

Total 1,959,651 records* 


*Valid records in original data sets obtained from National Space 
Science Data Center. Not all records contain both Ne and Te values. 




Local Time-Solar Zenith Angle 


Current database 



Local Time-Solar Zenith Angle 



o 50 100 150 

SZA (deg) 

All Data fAE— C.-D. DE-2. Explorer 221 


Local Time-Solar Zenith Angle {350 km - 450 kn 



0 50 100 150 

SZA {deg) 

All Data OE-C.-D. DE-2. Exolorer 221 
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Te, Ne vs Altitude 


ISS 

altitude 


Electron Temperature -Altitude Electron Density— Altitude 
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Ne, Te Scatter Plots 
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Te,Ne Scatter Plot 
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200 km < z < 1000 km 


~1.07xl0 6 values 

-90.0° < latitude < +90.0° 0° < SZA < 180° 




ISS Filtered Ne, Te Scatter Plots 
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Dole Ne (#/Vn3) 


IRI vs Data 


|lat| < 51.6 350 km < z < 450 km 



IRI 2001 Ne (§/ m3) 

Blue: mean IRI value in bin 
Green: +/- 68% Yellow: +/- 95% 
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International Reference Ionosphere (IRI) 


Environment variability: 

- 13 month FI 0.7 (solar EUV proxy) 

- Ap (2000) 

- Climatological model 

- Not intended to reproduce individual events 

Diurnal, spatial variations 
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• Good for design studies 

• Empirical model is fast 
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3 a Calculation 


- Bin satellite data based on increments of corresponding IRI model Ne values (logarithmic bins) 

• Within each bin, compute IRI Ne (or Te) in center of bin, IRI = (IRI2-IRIl)/2 

• Compute total number of data points in bin = m 

- Divide data in bin into two populations: 

• High Population: Data values > IRI Ne and Te within each bin 

• Low Population: Data values < IRI Ne and Te within each bin 

- Sequence all data point values (either Ne or Te) from lowest to highest = 1 to m 

- Define data (Ne or Te) points in bin which are below IRI (Ne or Te) 

• Number of points below IRI is ‘n’ and points range from 1 to ‘n’ 

• Point ‘n’ is highest data point (Ne or Te) which is less than IRI 

- Define data (Ne or Te) points in bin which are above IRI (Ne or Te) 

• n+l to m, where point n+1 has lowest value in upper bin above IRI 

- Computation of al, g2, g3 are based on number of points: 

• Compute upper bin oO, alu, a2u, o3u values (aO = 0) from values at percentiles of 0%, 68%, 95%, and 99.9% for 


upper bin 


au 


P = Point count from Bottom 


5Ne=Ne(P)-Ne(n+l) 


5Te=Te(P)-Te(n+l) 


0 

68% 

95% 

99.9% 


n+l 

0.68(m-n)+n 

0.95(m-n)+n 

0.999(m-n)+n 


SNe(P) 

SNe(P) 

8Ne(P) 

5Ne(P) 


8Te(P) = aOu’ 
8Te(P) = alu’ 
8Te(P) = a2u’ 
8Te(P) = g3u’ 


3 a Calculation 


- Computation of ol, cr2, a3 are based on number of points (continued): 

• Compute lower bin oO, Glu, g2u, g3u values (gO = 0) from values at percentiles of 0%, 68%, 95%, and 99.9% for 
lower bin 


ai 

P = Point count from Bottom 

8Ne=Ne(P)-Ne(n+l) 

6Te=Te(P)-Te(n+l) 

0% 

1.0(n) 

SNe(P) 

8Te(P) = aOP 

68% 

0.32(n) 

8Ne(P) 

8Te(P) = all’ 

95% 

0.05(n) 

SNe(P) 

8Te(P) = a2P 

99.9% 

0.001(n) 

8Ne(P) 

8Te(P) = a31’ 


- Final Table for each bin: 

m,n Number of points for Ne and Te 
IRI (Ne orTe) 
g values for Ne or Te 


aOu 

aOl 

alu 

all 

a2u 

a21 

a3u 

a31 


Variability is Ne or Te = IRI + a 


Dole Ne (#/rn3) 


IRI vs Data 


|lat| < 51.6 350 km < z < 450 km 




13 


Blue: mean IRI value in bin 
Green: +/- 68% Yellow: +/- 95% 


Algorithm 2, No AE— E 



Red: +/- 99.9% 
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Eclipse Exit Te 
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FI 07 Index Ne (#/m3) Te (eV) 


Solar Cycle Variations of Eclipse Exit Te 


Comparison ot Arecibo ISR Te.Ne and IRI Te,Ne (for ISS Charging Events) 
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Ne, Te Correlation 


Arecibo Te vs Ne 
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Current work 

- New data 

- Discussion 

• Documentation 

• 3 sigma 

• Schedule 



Old Data 


Satellite Data (in house) 


AE-C 

150 km 

X 

4300 km 

X 

68.1 deg 

AE-D 

154 km 

X 

13816 km 

X 

90.1 deg 

AE-E 

156 km 

X 

2983 km 

X 

90.1 deg 

DE 2 

300 km 

X 

1000 km 

X 

90 deg 



Incoherent Scatter Radar Data (in house) 
Arecibo (T/R) 18.3N -66.8E 85-99 
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New Data 


Satellite Data (in house) 


AE-C 

150 km 

X 

4300 km 

X 

68.1 deg 

AE-D 

154 km 

X 

13816 km 

X 

90.1 deg 

AE-E 

156 km 

X 

2983 km 

X 

90.1 deg 

DE 2 

300 km 

X 

1000 km 

X 

90 deg 

Explorer 22 

889 km 

X 

1081 km 

X 

79.7 deg 

Explorer 3 1 

500 km 

X 

3000 km 

X 

79.8 deg 

Aeros-A 

223 km 

X 

867 km 

X 

96.9 deg 

Alouette 2 

505 km 

X 

2987 km 

X 

79.8 deg 

Aureol 3 

380 km 

X 

1920 km 

X 

82.6 deg 

Hinotori 

548 km 

X 

603 km 

X 

3 1 .3 deg 

ISIS-1 

578 km 

X 

3526 km 

X 

88.42 deg 


Satellite Data (in process) 

DMSP 835 km x 835 km x 98.2 deg 

- 2000, 2001 , 2002 in hourse 

- 1998, 1999 ready for downloading 

CHAMP 454 km circular x 87 deg 

- AFRL is reprocessing using updated 
algorithms, available -August 2003 


Incoherent Scatter Radar Data (in house) 


Arecibo (T/R) 

18.3N 

-66.8E 

66,71-00 

St. Santin (T) 

44.6N 

2.2E 

66-86 

-Nancay (R) 

47.4N 

2.2E 


— Mende (R) 

44.5N 

3.45E 


— Monpazier (R) 

44.7N 

0.8E 


Millstone Hill (T/R) 

42.6 N 

-71.5E 

78-00 

Jicamarca (T/R) 

11.95N 

-76.0E 

66,86,99-00 

Irkutsk (T/R) 

52.85N 

103.1E 

99-02 

Karkov (T/R) 

50.0N 

36.2E 

02 

Chatanika (T/R) 

65. IN 

147.4W 

76-81 

Sondrestrpm 

67.0 N 

-51.0E 

83-02 

EISCAT,Troms0 (T/R) 

69.6 N 

19.2E 

84-99 

—Kiruna, Sweden (R) 

67.87N 

20.43E 


— Sodankyla, Finland (R) 

67.37N 

26.63E 


EISCAT, Svalbard (T/R) 

78N 

15.7E 

98 
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Altitude (km) 


FI 07 Index 


New Database 
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Altitude (km) 


Forward Work 


• Satellite data analysis 

- Complete acquisition of DMSP, CHAMP data 

- Complete processing satellite data holdings into study sets 

- Incorporate new data into statistics 

- Re-derive na statistics 

• Both original Ne, Te independent values 

• [Ne,Te] correlated statistics e.g., explore statistics of (Ne/Te) or similar parameters 

• Incoherent scatter radar data analysis 

- Data 

• Process from CEDAR Database CBF format to ascii (complete) 

• Process ascii to study sets (in process) 

• Modify analysis software to include radar data 

• Incorporate ISR data into statistics 

• Re-derive na statistics for ISR data 

- Other 

• Contact Pi’s for radar data (obligation of using CEDAR Database) 

• MU radar data? 


- CEDAR lists velocity data, but no NTV, <500 km. 

- Contact S. Fukao (PI). 




• Check CEDAR database for updated files (B. Emory/HAO indicated new data, May 2003) 
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Forward Work 


• Variability Model 

- Provide updated variability model to Boeing 

• Do we want options for statistics of (1) satellites, (2) ISR, and (3) both? 

- Coordinate model format with H. Barsamian for ease of incorporation into PIM 

• Documentation 

- Write report documenting data analysis, results 

• Analysis of satellite, radar data 

• Geophysical conditions that may yield high Ne, low Te (input from S. Koontz). 

- Peer review 

• Submit publications describing variability work to peer reviewed journal (with 
working group as co-authors) 

• Presentation at Spacecraft Charging conference? 

• Invitation to present variability work at 2003 International Reference Ionosphere 
Workshop (October 2003) 

- Papers published in Advances in Space Research 
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Overview 


• Ground based radar 

- Sites, database 

- Examples 

- Comparison with FPMU data 

• Low electron temperature discussion 


Review of Ground Based Radar Data 


ISR Sites 



•6 sites within ±51.6° latitude 

•NSF facilities at Millstone, Arecibo, Jicamarca best 
opportunity for our use 
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Operational ISR Sites 


Site 

Lat 

Lon 

Type 

Contact 

Jicamarca, Peru 

-11.95 

-76.0 

ISR 

Kudeki/Hysell 

Arecibo, Puerto Rico 

18.3 

-66.8 

ISR 

Q. Zhou 

Millstone Hill, Massachusetts 

42.6 

-71.5 

ISR 

Foster/Zhang 

Kharkov, Ukraine 

50.0 

36.2 

ISR 

V. Taran 

Irkutsk, Russia 

52.85 

103.1 

ISR 

G. Zherebtsov 

Sondrestrom, Greenland 

67.0 

-51.0 

ISR 

J. Thayer 

EISCAT Sodankyla, Finland 

67.4 

26.6 

ISR 

A. Van Eyken 

EISC AT Kiruna, Sweden 

67.9 

20.4 

ISR 

A. Van Eyken 

EISCAT Tromso, Norway 

69.6 

19.2 

ISR 

A. Van Eyken 

EISCAT Svalbard (ESR), Norway 

78.2 

15.7 

ISR 

A. van Eyken 

MU, Japan 

34.8 

136.1 

ISR 

S. Fukao 
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Millstone Hill ISR-ISS Geometry 
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7 



600 


11 Apr 22:42:39 UT 


(close approach time) 


2001/4/11 2 2 33; 36 - 2251;36 UT 



Lot Rude (de<g) 

(o) 


6001 i ■■■ i ■■ ■ 

500 ■ ■ 

^ 400 • 

-X 

* 300 ■ \ 

2 

% 200 ■ 

100 • 

o . . . . . . . i . > . 

-loo -aa —60 -40 -20 
LonqRude (deg) 

<b) 



8 


11 Apr 22:42:39 UT 


2001/4/11 2233:36 - 2251:36 UT 

1 n n e>r\n — 1 1 


io 9 io 1 c lo 1 ^ o 1 ^ a 1 3 

Ne (jf/m 3 ) 


eoa 

n ■ "■ ■'■■■1-- 

eoa ' 

^ 500 


50Q ; 

c 


L. 



' b — •* 

* 4oa 

* 

■g 400 ; 

□ 

m r 

□ 

< 300 

* . 

i • ' . t 

1 * • • 

< 30Q : 

20a 

1 * *• • 1 
i • /„ | 

200 . 


ioa 


1000 
Te (K) 



10000 


11 Apr 22:42:39 UT 


MILLS™ HILL RADAR 01/04/11 


Ne 



ISS intersects 
elevation scan at: 

42.8°N 

398.0 km 



25 30 35 40 45 50 55 60 65 


Lat i tude <deg) 



25 30 35 40 45 50 55 60 65 


L at i tude <deg) 



2000 4000 



CD 

CD 

in 


1 ' 

•o 

D CD 

■H © 
£ ” 


© 

© 



22:41:41-22:43:09 UT 

Ti 


25 30 35 40 45 50 55 60 65 


Lat i tude <deg) 



1500 2000 2500 3000 


Vi 


25 30 35 40 45 50 55 60 65 

Lat i -tude (deg) 

-600-400-:l?y^Sf 400 600 




12 Apr 03:29:00 UT 


2001/4/12 0323;42 - 0341-42 UT 

6QO 




(b) 


-100 -BO -6Q -40 -20 

Longrtude (deg) 




11 


12 Apr 03:29:00 UT 


£ 


T> 

□ 


600 

500 


< 300 


2001/4/12 032 3; 4 2 - 0341;42 UT 

■■■■ i J n n cr\n ■ — i 


1 — ■ ■"■“■I 1 

600 

. .*<• 

'£ 500 

-:r 

V- 


* 400 
□ 

* % 

< 300 

• . • i 

■! J ■ ...» 

200 


200 

10 9 10 1C lo 1 W^O 13 


Ne (tf/m 3 ) 


100 


1000 
Te (K) 


10000 


12 Apr 03:29:00 UT 


MILLSTONE KELL RADAR 01/04/12 

He 


03:17:49-03:19:33 UT 

Ti 



ISS intersects J 
elevation scan at: ■§ 

-f 


47.2°N 


399.1 km 


25 30 35 40 45 50 55 60 65 


L at i tude Cdeg) 



10.5 


ll’.0 


11.5 


12.0 



Lat i tude <deg) 



1000 2000 3000 




25 30 35 40 45 50 55 60 65 


Hi 

1000 


Lat i tude <deg> 



1500 2000 2500 




13 



Statistical Comparisons of FPP and ISR Results 


Encounter Date/Time 

AT isr : 18 min AT FPP : 

2001/4/11 2233:36- 2251:36 UT 
2001/4/12 0006:36 - 0024:36 UT 
2001/4/12 0143:48- 0201:48 UT 
2001/4/12 0323:42- 0341:42 UT 
2001/4/12 0500:00- 0518:00 UT 
2001/4/13 0043:48- 0101:48 UT 
2001/4/13 0224:00 - 0242:00 UT 
2001/4/13 0401:48 - 0419:48 UT 


FPP ISR 


Ne 

°Ne 

Ne 

^Ne 


(#/m 3 xlO 11 ) 


9.6 min 

x 2 < 0.05 

3.58 

5.09 

3.32 

3.93 

19.00 

4.69 

25.88 

142.21 

0.94 

0.29 

3.08 

8.10 

1.26 

0.98 

1.03 

1.21 

1.87 

1.47 
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Reanalysis of this data in process 
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Low Electron Temperature Discussion 
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Arecibo (18°N) Ne, Te 
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ISR Survey Data 
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ISR Survey Data 


St. Santin, France 
1966 - 1979 
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ISR Survey Data 


St. Santin, France 
1980 - 1986 





ISR Survey Data 
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ISR Survey Data 
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ISR Survey Data 
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ISR Survey Data 


Millstone Hill, MA 
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ISR Survey Data 


Millstone Hill, MA 
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Te.Ne Scatter Plot : 350km[v Altitude <450km 
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Forward Work 


• Process historical data sets from CEDAR Database for inclusion in 
variability study 

- Filter for ISS pre-sunrise conditions 

- Solar cycle pre-sunrise Te, Ne variation 

- Discuss ISR data analysis with Pi’s 

• Determine if any night ISR Te, Ne measurement will constrain ionosphere 
model sufficiently well to use for all eclipse exits 

- Compare available ISR world day data from multiple sites, satellite data 

• Re-analysis of FPP and Millstone Hill ISR data 

- L. Neergaard completed a more sophisticated model of Millstone radar beam 
including locations of all radar gates 

- ISS trajectory through Millstone Hill radar beam and identification of intersection 
of ISS trajectory with radar gates 

- Need to complete comparison of radar Ne, Te data from range gates that intersect 
ISS orbit with FPP Ne, Te values 
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Forward Work 


• Analysis of ionosonde Ne values with FPP overpasses 

- Useful to test process for comparing FPMU density values with ionosonde data 

- NOAA SPIDR database updated with AFRL ionosonde network data on 5-15 
minute time scale 

• Issue of h m F 2 retrieval accuracy for automated ionogram scaling used to fill 
SPIDR database (discussion with NOAA/SEC personnel at 2003 Space 
Weather Week). 

- Contact R. Vieriek/NOAA/SEC to determine what has been done on historical data 
analysis 

- AFRL to use ionosonde data in validation of DMSP satellite electron density 
retrieval from SSULI airglow limb observations 

• B. Reinisch/Lowell Observatory to manually check Air Force network 
ionograms, scaling 

• DMSP launch with SSULI is TBD. 

• ISS PCU-ISR radar correlation (S. Koontz) 

- Compare list of ISS PCU operation dates with ISR records to locate coincident 
observations during ISS overflights of radar sites 
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Forward Work 


• Write proposals for NSF ISR radar support of FPMU on-orbit check 

- Sell as science program validation of new charging model for high voltage solar arrays? 

- Millstone Hill 

• UHF Incoherent Scatter Radar requests submitted ideally 6 weeks before the desired experiment 
date 

• Coordinated observations between UHF Radar and other instruments submitted ideally 3 months 
before the desired experiment date. 

- Arecibo 

• Proposals due 2 June 2003 for observations from 1 October 2003 - 31 May 2004 

• Need to submit now if FPMU is flown to ISS in March 2004 (this weekend?!) 

- Jicamarca 

• Scheduling experiments at Jicamarca is still handled in an informal way. Contact D. Farley at 
Cornell and/or Ronald Woodman at JRO for details. 

• If you plan to run at night or during Friday-Sunday, you should be prepared to pay overtime charges 
to the staff members. These charges might add up to $30-50/hour, depending upon the number of 
people involved. Staff members are generally very happy to work overtime because payments 
represent a substantial boost to their income. 

• Cost (NSF supported research): no charge for observing time, other than overtime for staff. 

• Cost(non-NSF supported programs): 

- There is a charge of approximately $8000 per day of observing for isolated experiments 

- For longer, on-going programs not supported by the NSF, special arrangements must be made with the 
Observatory. 



